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Abstract 

Heavy ion collisions at RHIC and at the LHC can cre- 
ate the quark-gluon plasma, a state of matter at very high 
temperatures. Among a plethora of particles that are pro- 
duced in these collisions, also light is emitted throughout 
the evolution of the plasma. 

In this talk, the properties of this light are discussed and 
related to recent efforts towards shorter and more energetic 
photon pulses in laser physics. In particular, the time evo- 
lution of high-energy photons is studied. These photons 
originate from Compton scattering of gluons and quark- 
antiquark annihilation in the plasma. Due to the internal 
dynamics of the plasma, double pulses at the yoctosecond 
time scale can be generated under certain conditions. Such 
double pulses may be utilized for novel pump-probe exper- 
iments at nuclear time scales. 

INTRODUCTION 

The year 2010 marks the fiftieth anniversary of the in- 
vention of the laser Since its invention, not only did the 
intensity steadily increase, but also the pulse duration be- 
came shorter and shorter In fact, pulse duration and in- 
tensity of lasers (or derived coherent radiation bursts) turn 
out to be correlated over a large range of energies and time 
scales m. This observation provides a good motivation to 
present at a conference on Physics in Intense Fields (PIF 
2010) a study of the shortest possible light flashes that can 
be produced in experiments, which is heavy ion collisions 
that produce the quark-gluon plasma (QGP) |2 1. 

A few selected milestones of the development of short 
laser pulses are depicted in Fig. [1] Each of the time scales 
has enabled to access new systems: At the picosecond 
to femtosecond timescale, femtochemistry allows for the 
time-resolved study of chemical reactions [3 |. For pump- 
probe spectroscopy, it is essential to have two short pulses 
in close succession: The first laser pulse triggers a chem- 
ical reaction, which may involve an excited state and var- 
ious short-lived transitions, while the second pulse takes a 
snapshot of the intermediate state. By varying the interval 
between the two pulses, the time-evolution of a chemical 
reaction can be studied. 

At shorter timescales, attosecond science is the window 
to capturing electron motion in molecules and atoms |4|. 
High-order harmonics of femtosecond laser radiation have 
been shown to be sources of trains of attosecond extreme- 
ultraviolet pulses |5J that can be used to produce single 
attosecond soft X-ray pulses |i6J. For example, such sin- 
gle attosecond X-ray bursts Q have applications in molec- 
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Figure 1: History of laser pulse duration. Marked are 
selected milestones of technologies that allowed to de- 
crease the laser pulse duration, like Neodymium glass laser 
(Nd:glass), Continuous Wave Dye laser (CW Dye), Col- 
liding Pulse-Mode locked dye laser (CPM), Titanium sap- 
phire laser (Ti:sapphire), or High-Harmonic Generation 
(HHG). For comparison, also the lifetime of the Quark 
Gluon Plasma (QGP) is indicated as it is produced nowa- 
days in heavy ion colliders like RHIC or LHC. 

ular imaging fSj, quantum control ||9|, or Raman spec- 
troscopy ifTol . By introducing a controlled delay between 
two such peaks, the dynamics of electron systems could be 
studied using pump-probe techniques |6|. Such techniques 
also allow for the direct time resolution of many-body dy- 
namics, like the observation of the dressing process of 
charged particles |fTT|. In the zeptosecond regime, nuclear 
processes may become accessible [12|. It has been sug- 
gested that zeptosecond pulses could be created via nonlin- 
ear Thomson backscattering 1 13 , 14 1, or by employing rela- 
tivistic laser-plasma interactions lfT5]iT6l . A possible detec- 
tion scheme for the characterization of short 7-ray pulses 
of MeV to GeV energy photons down to the zeptosecond 
scale has been proposed in Ref. IfTTI . 

At even shorter timescales, double pulses of yoctosecond 
duration of GeV photon energy could be created from the 
quark-gluon plasma in non-central heavy ion collisions ||2j| . 
It turns out that the emission envelope can depend strongly 
on the internal dynamics of the QGP. Under certain con- 
ditions, a double peak structure in the emission envelope 
may be observed. This could be the first source for pump- 
probe experiments at the yoctosecond timescale. The de- 
lay between the peaks is directly related to the isotropiza- 
tion time, and the relative height between the peaks can be 
shaped by varying photon energy and emission angle. Such 
pulses could be utilized, for example, to resolve dynam- 
ics on the nuclear timescale such as that of baryon reso- 
nances ifTSl . Conversely, a time-resolved study of the emit- 



ted photons could provide a window to the internal QGP 
dynamics throughout its expansion. 

QGP PHOTON PRODUCTION 

The Large Hadron ColUder (LHC) at CERN has be- 
gun its heavy-ion program in November 2010 ||T9l . just a 
few weeks before the PIF 2010 conference. With center- 
of-mass energies of 2.76 TeV per nucleon, the collisions 
of lead ions produce hotter and denser plasmas than pre- 
viously achievable at the Relativistic Heavy Ion Collider 
(RHIC), where gold ions were used. The temperatures 
reached in these collisions are so high that the constituents 
of atomic nuclei, the neutrons and protons, are split into 
their constituents, the quarks and gluons. The interest in 
the QGP stems not least from the fact that it is believed to 
have filled the entire universe during the first few microsec- 
onds after the Big Bang. 

In heavy-ion collisions, the QGP is produced up to the 
size of a nucleus (^15 fm) for a duration of a few tens of 
yoctoseconds (lys = 10^^"' s). The plasma is produced 
initially in a very anisotropic state, and reaches a hydro- 
dynamic evolution through internal interactions only after 
some isotropization time riso (see Fig.|2l). Among the many 
particles that are produced, also photons of a few GeV en- 
ergy are emitted ll20l |2T1 . 

It was one of the early surprises that the observed parti- 
cle spectra turned out to agree well with ideal hydrodynam- 
ical model predictions |22|. This led to the initial assump- 
tion that isotropization times may be as low as Tjso ~ 1 ys. 
However, it has been pointed out in the mean-time that vis- 
cous hydrodynamic models are still consistent with RHIC 
data if isotropization times as large as Tiso w 7 ys are as- 
sumed li23l . even if the expansion before isotropization 
is assumed to be collisionless ("free streaming"). This 
should be compared to the lifetime of the QGP, which could 
amount to 15 ys at RHIC, and which could be as large as 
25 ys at LHC. 

Figure 12 a-c) shows a schematic view of the collision of 
two heavy ions. The ions are illustrated as relativistically 
contracted pancakes. In general, two ions will not collide 
head-on-head, but will be displaced by an impact parame- 
ter b. Direct photons are emitted from the expanding QGP 
throughout its lifetime |24|. The energy spectrum of the 
emitted photons extends to the GeV range, and the upper 
limit for the temporal duration of the GeV photon pulse is 
determined by the expansion dynamics of the QGP, which 
leads to yoctosecond pulses. At the initial stage of the colli- 
sion, even before the plasma is created, prompt photons are 
emitted from nucleon-nucleon collisions in all directions, 
see Figs. 12 a) and (d). For an intermediate time after the 
collision shown in Fig|2b), a momentum anisotropy oc- 
curs due to the longitudinal expansion of the plasma: Those 
particles that originally had momentum components in for- 
ward or backward direction along the beam axis leave the 
central region quickly, so that mainly particles with trans- 
verse momenta remain in the plasma. High-energy pho- 




Figure 2: Early stages of a high-energy collision, involv- 
ing pre-equilibrium (first two columns) and equilibrated 
QGP phases (last column). Parts (a)-(c) show three snap- 
shots in time in position space. Shown are the two rela- 
tivistically contracted colliding ions that create the quark- 
gluon plasma in the overlap region. Curly arrows de- 
note photon emission and semicircles the detectors. Parts 
(d)-(f) are corresponding pictorial representations of the 
plasma in momentum space. In an intermediate stage of 
the pre-equilibrium phase, the momentum distribution is 
anisotropic, resulting in a change in the angular photon 
emission pattern that can give rise to double-peaked pho- 
ton pulses 12) ■ 

tons that are created in the plasma through Compton scat- 
tering or quark-antiquark annihilation carry preferentially 
the momentum of the original participants of the collision. 
Because of the momentum anisotropy of the quarks and 
gluons within the plasma, also the emitted photons are pref- 
erentially emitted perpendicular to the beam axis z, as in- 
dicated in Fig. |2te). Finally, in Fig. |2c), the system had 
time to isotropize due to collisions within the plasma. The 
photons will be emitted again in all directions, as shown 
in Fig. I2f), particularly also into the direction in which 
the photon emission was suppressed during the anisotropic 
stage. 

A photon detector placed towards the beam axis would 
therefore measure a time-dependent photon flux. Ideally, in 
the intermediate stage in Fig. I2e), the photon emission in 
the direction of the detector is suppressed highly enough, so 
that the photons emitted in the stages Fig.|2d) and Fig.|2f) 
are distinct enough in time to form two separate photon 
pulses. In order to quantitatively describe the pulse enve- 
lope of the emitted photons, detailed calculations are nec- 
essary. Such a calculation has been performed in Ref. 
which was based on a one-dimensional expansion model 
by Bjorken ll25l . This model assumes boost-invariant evo- 
lution of the quark-gluon plasma in a central region of the 
collision. The detector is placed away from the beam axis 
by an angle 6 within the reaction plane. Direct photons in 
the GeV energy range are emitted from the pre-equilibrium 
and equilibrated phases of the QGP. The leading contribu- 
tion to the photon production rate R originates from quark- 
gluon Compton scattering and quark-antiquark annihila- 
tion. In principle, higher order soft scattering processes like 



bremsstrahlung or inelastic pair annihilation would have to 
be included as well, but their contributions become dom- 
inant only at lower energies, and are less important at the 
higher energies considered ll26l . 

For anisotropic momentum distributions, the photon pro- 
duction rate R has to be calculated numerically 1 2 1 ) . It 
depends on the temperature T of the medium, the pho- 
ton energy E and momentum k, the fine structure con- 
stant a, and the corresponding quantity for the strong force 
as (with h = c = = 1). The rate further de- 
pends on the anisotropy, which is described by a parameter 
C = (pt) / (2 {pl}) ~ 1 that relates the mean longitudi- 
nal and transverse momenta and ll26l l27l . To inte- 
grate this rate over time, a time evolution model for the pre- 
equilibrium and equilibrated QGP has been used [27J. This 
model specifies the time evolution for the energy density 
E = S(t), for the hard scale phard — Phard(''') (which cor- 
responds to T in the isotropic case), and for the anisotropy 
parameter ^ = ^(r) as a function of the proper time r. 
Qualitatively, the model follows the evolution as outlined 
in Fig. 12] For early times, a free streaming phase lets the 
anisotropy grow. At late times, the system converges to 
an ideal hydrodynamic phase with vanishing anisotropy. 
These two phases are linked by a smooth transition which 
is controlled by additional model parameters. Thermahza- 
tion and isotropization happen concurrently in this model, 
Tthorm — T\^o- The model is thus able to cover both, the 
pre-equilibrium phase and the equilibrated QGP phase of 
the expanding plasma. 
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Figure 3: Photon emission rate as a function of detector 
time r. Solid blue lines show a large isotropization time 
Tiso = 6.7 ys while dashed red lines correspond to a short 
isotropization time Tiso = 0.3 ys. Parts (a) and (b) dis- 
play emission at midrapidity {B = 7r/2) for a central col- 
lision with impact parameter 6 = 0. Parts (c)-(f) show 
double-peaked photon pulses obtained for b = 9.2 fm or 
12.2 fm, and the vertical dotted line indicates the position 
of the larger Tjso — 6.7 ys. In parts (c) and (d), the detector 
direction is 6* = 7r/4, in (e) and (f), it is 6* = 7r/8 10. 



The numerical parameters suitable for calculating LHC 
parameters are given as follows: The initial temperature 
is assumed to be Tq = 845 MeV with a formation time 
To = 0.3 ys. The critical temperature, where the QGP 
ceases to exist, is taken as Tc = 160 MeV. The isotropiza- 
tion time is varied in the range Tiso — to t^so — 6.7 ys, 
assuming free-streaming at early times. Both possibilities 
are not yet ruled out by RHIC data. In order to ensure 
fixed final multiplicity, the initial conditions are adjusted as 
a function of Tiso such that the same entropy is generated 

as for Tiso = Tq. 

For central collisions with emission angle orthogonal to 
the beam axis (9 = tt/2), a typical time evolution of the 
photon emission rate is depicted in Figs.[3la) and[3lb). At 2 
GeV energy, the photon production from the QGP at midra- 
pidity is 3 to 4 times as large as the production from the 
initial collisions. It is roughly 6 times as large as the pro- 
duction from the hadron gas [24|. At 3 GeV energy, QGP 
photon production is even 50 times larger than the produc- 
tion from the hadron gas. In the Figs. |3] the origin of the 
abscissa is the time when a photon emitted from the cen- 
ter of the collision arrives at the detector Photons arriving 
earlier originate from a part of the QGP that is closer to the 
detector. The pulse shape is mainly determined by the ge- 
ometry of the lead ion with radius 7.1 fm. Any structure on 
the yoctosecond timescale is blurred simply by the time for 
fight to traverse the QGP. 

This limit can be overcome in the following ways: By 
considering non-central collisions with impact parameter 
b, the physical extent of the QGP is reduced. Also, an opti- 
mization of the detection angle can minimize the traveling 
time through the plasma. In forward direction, the initial 
shape of the QGP is Lorentz-contracted, and light leaves 
this initial region quickly. This is partially spoiled due to 
the QGP expansion in the same direction. Thus intermedi- 
ate emission angles are most promising for which the QGP 
appears partly Lorentz contracted but does not expand to- 
wards the detector 

Figures [3jc)-13f) show the photon emission in the direc- 
tions 9 = tt/A and 6 — tt/8. For large impact parameters 
b = 9.2 fm and b — 12.2 fm, a striking double-peak struc- 
ture appears. The minimum between the two peaks corre- 
sponds roughly to maximum anisotropy within the plasma. 
This follows from the fact that the photon emission rate 
is suppressed for larger values of f and smaller values of 
9 Il26l . The distance between the two peaks is approxi- 
mately governed by the isotropization time Tiso, indicated 
by the dotted line in Figs. [3jc)-[3]!f)- The first peak corre- 
sponds to photons emitted from the blue-shifted approach- 
ing part of the QGP, while the second peak corresponds to 
a slightly red-shifted and time-dilated receding tail of the 
plasma. For a short isotropization time Tiso = tq (dashed 
lines) the separation into two peaks does not occur There- 
fore, this effect depends delicately on the QGP dynamics. 

There are a couple of caveats to this model calculation: 
Apart from the photons originating from the QGP, in an 
actual experiment there is a background of photons from 



various sources f2S\. These include photons produced by a 
jet passing through the QGP 1.20 j , and could dominate the 
effect that is expected from the QGP alone. Since these 
photons are produced on a similar yoctosecond timescale, 
they would modify the pulse shape on this timescale. Pho- 
tons produced from the initial collisions of the two ions 
can be of comparable size in intensity, but they would 
only enhance the first peak of the double peaks depicted 
in Fig. [3lc)-|3jf). Other background photons are produced 
in the decay of pions at later stages of the collision. Since 
these are produced at much later time scales, they would 
not modify a time structure on the yoctosecond timescale. 
It would be necessary to take these various photon sources 
into account in order to obtain a quantitative prediction of 
the expected pulse structure. 

What kind of properties should the detector have in or- 
der to resolve the double pulses? The time-integrated high- 
energetic photons are already routinely detected in experi- 
ments |29|. In the GeV energy range, the photon produc- 
tion rate is of the order of a few photons per collision [28 1. 
Note that a single GeV photon pulse of 10 ys duration cor- 
responds to a pulse energy of only about 100 pj, but to a 
power of 10 TW. The photon yield could be enhanced by 
considering lower energy photons, but this would also in- 
crease the number of unwanted background photons. Alter- 
natively, the photon yield could be enhanced by increasing 
the collision energy. This may have the additional benefit 
of increasing the relative importance of the contribution of 
thermal photons compared to other kinds of photons [,28 1 . 
The emission envelope is influenced by the geometry, emis- 
sion angle, and internal dynamics like the isotropization 
time of the expanding QGP. The double-peak structure de- 
scribed may emerge in non-central collisions at an emis- 
sion angle close to forward direction, assuming that the 
isotropization time is large. In order to detect such short 
pulses, new detection schemes would be required. Exist- 
ing tools and ideas from attosecond metrology, like pump- 
probe experiments or spectroscopy techniques, may turn 
out to be appropriate candidates to be scaled to zepto- or 
yoctosecond duration [T TTl. Also, an experimental deter- 
mination of the photon emission envelope would serve as 
an additional probe of the internal dynamics of the QGP, 
for example by measuring its isotropization time. 

CONCLUSIONS 

The steady progress in laser physics over the last decades 
towards higher intensities and shorter pulse lengths pro- 
vides strong motivation to study systems that produce ex- 
tremely short flashes of light. The QGP is an ideal candi- 
date because it exhibits non-trivial dynamics on the yoc- 
tosecond timescale, during which GeV photons are emit- 
ted. Under certain conditions, a double-peak structure may 
be produced. This could eventually lead to novel pump- 
probe experiments at the GeV energy scale. Alternatively, 
measuring the temporal shape of the photon emission enve- 
lope, dynamic properties of the QGP, like its isotropization 



time, could be probed experimentally. 

ACKNOWLEDGMENT 

I would like to thank my collaborators J. Evers, 
K. Z. Hatsagortsyan, and C. H. Keitel for guidance and 
fruitful discussions that led to the work that has been pre- 
sented |l2][l7l. I would further like to thank the organizers 
of flie PIF 2010 for flieir kind hospitality. 

REFERENCES 

[1] G. Mourou and T. Tajima. Science, 331, 41 (2011). 

[2] A. Ipp, C. H. Keitel, and J. Evers. Phys. Rev. Lett., 103, 
152301 (2009). 

[3] A. H. Zewail. Femtochemistry: Ultmfast Dynamics of the 
Chemical Bond, World Scientific (1994). 

[4] M. Drescher, et al. Science, 291, 1923 (2001). 

[5] P. Paul et al. Science, 292, 1689 (2001). 

[6] Y. Silberberg. Nature, 414, 494 (2001). 

[7] T. Baeva, S. Gordienko, and A. Pukhov. Laser Part. Beams, 
25, 339 (2007). 

[8] M. Lein. Physical review letters, 94, 053004 (2005). 

[9] H. Rabitz, et al. Science, 288, 824 (2000). 

[10] N. Dudovich, D. Oron, and Y. Silberberg. Nature, 418, 512 
(2002). 

[11] R. Huber, et al. Nature, 414, 286 (2001). 

[12] C. Golabek and C. Simenel. Phys. Rev. Lett., 103, 042701 
(2009). 

[13] P Lan et al. Phys. Rev. E, 72, 066501 (2005). 

[14] D. Kim etal. New J. Phys., 11, 063050 (2009). 

[15] S. V. Bulanov, T. Zh. Esirkepov, and T. Tajima. Phys. Rev. 
Lett., 91, 085001 (2003). 

[16] Y. Nomura etal. Nature Phys., 5, 124 (2009). 

[17] A. Ipp, ef a/. (2010). I arXiv: 1008.03551 . 

[18] M. Dugger et al. Phys. Rev., C76, 02521 1 (2007). 

[19] K. Aamodt et al. (2010). 

[20] R. J. Fries, B. MuUer, and D. K. Srivastava. Phys. Rev. Lett., 
90, 132301 (2003). 

[21] A. Ipp, etal. Phys. Lett., B666, 315 (2008). 

[22] P Huovinen et al. Phys. Lett., B503, 58 (2001). 

[23] M. Luzum and P Romatschke. Phys. Rev., C78, 034915 
(2008). 

[24] S. Turbide, R. Rapp, and C. Gale. Phys. Rev., C69, 014903 
(2004). 

[25] J. D. Bjorken. Phys. Rev, D27, 140 (1983). 

[26] B. Schenke and M. Strickland. Phys. Rev, D76, 025023 
(2007). 

[27] M. Martinez and M. Strickland. Phys. Rev. Lett., 100, 
102301 (2008). 

[28] S. Turbide, etal. Phys. Rev, C72, 014906 (2005). 

[29] S. S. Adler et al. Phys. Rev Lett., 94, 232301 (2005). 



